Signatures of Mutational Processes and Associated Risk Factors in Esophageal Squamous Cell Carcinoma: A Geographically Independent Stratification Strategy?

Cancer Research UKThis is the accepted manuscript. The final version is available at http://www.sciencedirect.com/science/article/pii/S001650851600370X

Visualization and analysis strategies for dynamic gene-phenotype relationships and their biological interpretation

The complexity of biological systems is one of their most fascinating and, at the same time, most cryptic aspects. Despite the progress of technology that has enabled measuring biological parameters at deeper levels of detail in time and space, the ability to decipher meaning from these large amounts of heterogeneous data is limited. In order to address this challenge, both analysis and visualization strategies need to be adapted to handle this complexity. At system-wide level, we are still limited in our ability to infer genetic and environmental causes of disease, or consistently compare and link phenotypes. Moreover, despite the increasing availability of time-resolved experiments, the temporal context is often lost. In my thesis, I explored a series of analysis and visualization strategies to compare and connect dynamic phenotypic outcomes of cellular perturbations in a genetic and network context. More specifically, in the first part of my thesis, I focused on the cell cycle as one of the best examples of a complex, highly dynamic process. I applied analysis and data integration methods to investigate phenotypes derived from cell division failure. I examined how such phenotypes may arise as a result of perturbations in the underlying network. To this purpose, I investigated the role of short structural elements at binding interfaces of proteins, called linear motifs, in shaping the cell division network. I assessed their association to different phenotypes, in the context of local perturbations and of disease. This analysis enabled a more detailed understanding of the regulatory mechanisms beyond the malfunctioning of cell division processes, but the ability to compare phenotypes and track their evolution was limited. Exploring large-scale, time-resolved phenotypic screens is still a bottleneck, especially in the visualization area. To help address this question, in the subsequent parts of the thesis I proposed novel visualization approaches that would leverage pattern discovery in such heterogeneous, dynamic datasets and enable the generation of new hypotheses. First, I extended an existing visualization tool, Arena3D, to enable the comparison of phenotypes in a genetic and network context. I used this tool to continue the exploration of phenotype-wide differences between outcomes of gene function suppression within mitosis. I also applied it to an investigation of systemic changes in the network of embryonic stem cell fate determinants upon downregulation of the pluripotency factor Nanog. Second, time-resolved tracking of phenotypes opens up new possibilities in exploring how genetic and phenotypic connections evolve through time, an aspect that is largely missing in the visualization area. I developed a novel visualization approach that uses 2D/3D projections to enable the discovery of genetic determinants linking phenotypes through time. I used the resulting tool, PhenoTimer, to investigate the patterns of transitions between phenotypes in cell populations upon perturbation of cell division and the timing of cancer-relevant transcriptional events. I showed the potential of discovering drug synergistic effects by visual mapping of similarities in their mechanisms of action. Overall, these approaches help clarify aspects of the consequences of cell division failure and provide general visualization frameworks that should be of interest to the wider scientific community, for use in the analysis of multidimensional phenotypic screens

Visualizing time-related data in biology, a review

Time is of the essence, also in biology. Monitoring disease progression or timing developmental defects are key aspects in the process of drug discovery and therapy trial. Furthermore, before deciphering the course of evolution of these complex processes, we need an understanding of the basic dynamics of biological phenomena that are often strictly time-regulated (e.g. circadian rhythms). With the advances in technologies able to measure timing effects and dynamics of regulatory aspects, visualization and analysis tools try to keep up the pace with the new challenge. Beyond the classical timeline plots, notable attempts at more involved temporal interpretation have been made in the recent years, but awareness of the available resources is still limited within the scientific community. Here we review some of the advances in biological visualization of time-driven processes and look at how they allow analyzing data now and in the future

Multi-scale characterisation of homologous recombination deficiency in breast cancer

BACKGROUND: Homologous recombination is a robust, broadly error-free mechanism of double-strand break repair, and deficiencies lead to PARP inhibitor sensitivity. Patients displaying homologous recombination deficiency can be identified using 'mutational signatures'. However, these patterns are difficult to reliably infer from exome sequencing. Additionally, as mutational signatures are a historical record of mutagenic processes, this limits their utility in describing the current status of a tumour. METHODS: We apply two methods for characterising homologous recombination deficiency in breast cancer to explore the features and heterogeneity associated with this phenotype. We develop a likelihood-based method which leverages small insertions and deletions for high-confidence classification of homologous recombination deficiency for exome-sequenced breast cancers. We then use multinomial elastic net regression modelling to develop a transcriptional signature of heterogeneous homologous recombination deficiency. This signature is then applied to single-cell RNA-sequenced breast cancer cohorts enabling analysis of homologous recombination deficiency heterogeneity and differential patterns of tumour microenvironment interactivity. RESULTS: We demonstrate that the inclusion of indel events, even at low levels, improves homologous recombination deficiency classification. Whilst BRCA-positive homologous recombination deficient samples display strong similarities to those harbouring BRCA1/2 defects, they appear to deviate in microenvironmental features such as hypoxic signalling. We then present a 228-gene transcriptional signature which simultaneously characterises homologous recombination deficiency and BRCA1/2-defect status, and is associated with PARP inhibitor response. Finally, we show that this signature is applicable to single-cell transcriptomics data and predict that these cells present a distinct milieu of interactions with their microenvironment compared to their homologous recombination proficient counterparts, typified by a decreased cancer cell response to TNFα signalling. CONCLUSIONS: We apply multi-scale approaches to characterise homologous recombination deficiency in breast cancer through the development of mutational and transcriptional signatures. We demonstrate how indels can improve homologous recombination deficiency classification in exome-sequenced breast cancers. Additionally, we demonstrate the heterogeneity of homologous recombination deficiency, especially in relation to BRCA1/2-defect status, and show that indications of this feature can be captured at a single-cell level, enabling further investigations into interactions between DNA repair deficient cells and their tumour microenvironment

Genomic and microenvironmental heterogeneity shaping epithelial-to-mesenchymal trajectories in cancer

The epithelial to mesenchymal transition (EMT) is a key cellular process underlying cancer progression, with multiple intermediate states whose molecular hallmarks remain poorly characterised. To fill this gap, we present a method to robustly evaluate EMT transformation in individual tumours based on transcriptomic signals. We apply this approach to explore EMT trajectories in 7180 tumours of epithelial origin and identify three macro-states with prognostic and therapeutic value, attributable to epithelial, hybrid E/M and mesenchymal phenotypes. We show that the hybrid state is relatively stable and linked with increased aneuploidy. We further employ spatial transcriptomics and single cell datasets to explore the spatial heterogeneity of EMT transformation and distinct interaction patterns with cytotoxic, NK cells and fibroblasts in the tumour microenvironment. Additionally, we provide a catalogue of genomic events underlying distinct evolutionary constraints on EMT transformation. This study sheds light on the aetiology of distinct stages along the EMT trajectory, and highlights broader genomic and environmental hallmarks shaping the mesenchymal transformation of primary tumours

Molecular effects of Lapatinib in the treatment of HER2 overexpressing oesophago-gastric adenocarcinoma.

BACKGROUND: Lapatinib, a dual EGFR and HER2 inhibitor has shown disappointing results in clinical trials of metastatic oesophago-gastric adenocarcinomas (OGAs), and in vitro studies suggest that MET, IGFR, and HER3 confer resistance. This trial applied Lapatinib in the curative neoadjuvant setting and investigated the feasibility and utility of additional endoscopy and biopsy for assessment of resistance mechanisms ex vivo and in vivo. METHODS: Patients with HER2 overexpressing OGA were treated for 10 days with Lapatinib monotherapy, and then in combination with three cycles of Oxaliplatin and Capecitabine before surgery. Endoscopic samples were taken for molecular analysis at: baseline including for ex vivo culture +/- Lapatinib to predict in vivo response, post-Lapatinib monotherapy and at surgery. Immunohistochemistry (IHC) and proteomic analysis was performed to assess cell kinetics and signalling activity. RESULTS: The trial closed early (n=10) due to an anastomotic leak in two patients for which a causative effect of Lapatinib could not be excluded. The reduction in Phosphorylated-HER2 (P-HER2) and P-EGFR in the ex vivo-treated biopsy demonstrated good correlation with the in vivo response at day 10. Proteomic analysis pre and post-Lapatinib demonstrated target inhibition (P-ERBB2, P-EGFR, P-PI3K, P-AKT, and P-ERK) that persisted until surgery. There was also significant correlation between the activation of MET with the level of P-Erk (P=0.0005) and P-PI3K : T-PI3K (total PI3K) ratio (P=0.0037). There was no significant correlation between the activation status of IGFR and HER3 with downstream signalling molecules. CONCLUSIONS: Additional endoscopy and biopsy sampling for multiple biomarker endpoints was feasible and confirmed in vitro data that MET is likely to be a significant mechanism of Lapatinib resistance in vivo.This research was funded by the Medical Research Council [Grant SK002].This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/bjc.2015.34

Whole-genome sequencing of nine esophageal adenocarcinoma cell lines.

Esophageal adenocarcinoma (EAC) is highly mutated and molecularly heterogeneous. The number of cell lines available for study is limited and their genome has been only partially characterized. The availability of an accurate annotation of their mutational landscape is crucial for accurate experimental design and correct interpretation of genotype-phenotype findings. We performed high coverage, paired end whole genome sequencing on eight EAC cell lines-ESO26, ESO51, FLO-1, JH-EsoAd1, OACM5.1 C, OACP4 C, OE33, SK-GT-4-all verified against original patient material, and one esophageal high grade dysplasia cell line, CP-D. We have made available the aligned sequence data and report single nucleotide variants (SNVs), small insertions and deletions (indels), and copy number alterations, identified by comparison with the human reference genome and known single nucleotide polymorphisms (SNPs). We compare these putative mutations to mutations found in primary tissue EAC samples, to inform the use of these cell lines as a model of EAC.This work was funded by an MRC Programme Grant to R.C.F. and a Cancer Research UK grant to PAWE. The pipeline for mutation calling is funded by Cancer Research UK as part of the International Cancer Genome Consortium. G.C. is a National Institute for Health Research Lecturer as part of a NIHR professorship grant to R.C.F. AGL is supported by a Cancer Research UK programme grant (C14303/A20406) to Simon Tavaré and the European Commission through the Horizon 2020 project SOUND (Grant Agreement no. 633974)

Immune Cell Abundance and T-cell Receptor Landscapes Suggest New Patient Stratification Strategies in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a molecularly and spatially heterogeneous disease frequently characterized by impairment of immunosurveillance mechanisms. Despite recent success with immunotherapy treatment, disease progression still occurs quickly after treatment in the majority of cases, suggesting the need to improve patient selection strategies. In the quest for biomarkers that may help inform response to checkpoint blockade, we characterized the tumor microenvironment (TME) of 162 HNSCC primary tumors of diverse etiologic and spatial origin, through gene expression and IHC profiling of relevant immune proteins, T-cell receptor (TCR) repertoire analysis, and whole-exome sequencing. We identified five HNSCC TME categories based on immune/stromal composition: (i) cytotoxic, (ii) plasma cell rich, (iii) dendritic cell rich, (iv) macrophage rich, and (v) immune-excluded. Remarkably, the cytotoxic and plasma cell rich subgroups exhibited a phenotype similar to tertiary lymphoid structures (TLS), which have been previously linked to immunotherapy response. We also found an increased richness of the TCR repertoire in these two subgroups and in never smokers. Mutational patterns evidencing APOBEC activity were enriched in the plasma cell high subgroup. Furthermore, specific signal propagation patterns within the Ras/ERK and PI3K/AKT pathways associated with distinct immune phenotypes. While traditionally CD8/CD3 T-cell infiltration and immune checkpoint expression (e.g., PD-L1) have been used in the patient selection process for checkpoint blockade treatment, we suggest that additional biomarkers, such as TCR productive clonality, smoking history, and TLS index, may have the ability to pull out potential responders to benefit from immunotherapeutic agents. // Significance: Here we present our findings on the genomic and immune landscape of primary disease in a cohort of 162 patients with HNSCC, benefitting from detailed molecular and clinical characterization. By employing whole-exome sequencing and gene expression analysis of relevant immune markers, TCR profiling, and staining of relevant proteins involved in immune response, we highlight how distinct etiologies, cell intrinsic, and environmental factors combine to shape the landscape of HNSCC primary disease

Mutational signature dynamics shaping the evolution of oesophageal adenocarcinoma

A variety of mutational processes drive cancer development, but their dynamics across the entire disease spectrum from pre-cancerous to advanced neoplasia are poorly understood. We explore the mutagenic processes shaping oesophageal adenocarcinoma tumorigenesis in 997 instances comprising distinct stages of this malignancy, from Barrett Oesophagus to primary tumours and advanced metastatic disease. The mutational landscape is dominated by the C[T > C/G]T substitution enriched signatures SBS17a/b, which are linked with TP53 mutations, increased proliferation, genomic instability and disease progression. The APOBEC mutagenesis signature is a weak but persistent signal amplified in primary tumours. We also identify prevalent alterations in DNA damage repair pathways, with homologous recombination, base and nucleotide excision repair and translesion synthesis mutated in up to 50% of the cohort, and surprisingly uncoupled from transcriptional activity. Among these, the presence of base excision repair deficiencies show remarkably poor prognosis in the cohort. In this work, we provide insights on the mutational aetiology and changes enabling the transition from pre-neoplastic to advanced oesophageal adenocarcinoma

Genomic hallmarks and therapeutic implications of G0 cell cycle arrest in cancer

BACKGROUND: Therapy resistance in cancer is often driven by a subpopulation of cells that are temporarily arrested in a non-proliferative G0 state, which is difficult to capture and whose mutational drivers remain largely unknown. RESULTS: We develop methodology to robustly identify this state from transcriptomic signals and characterise its prevalence and genomic constraints in solid primary tumours. We show that G0 arrest preferentially emerges in the context of more stable, less mutated genomes which maintain TP53 integrity and lack the hallmarks of DNA damage repair deficiency, while presenting increased APOBEC mutagenesis. We employ machine learning to uncover novel genomic dependencies of this process and validate the role of the centrosomal gene CEP89 as a modulator of proliferation and G0 arrest capacity. Lastly, we demonstrate that G0 arrest underlies unfavourable responses to various therapies exploiting cell cycle, kinase signalling and epigenetic mechanisms in single-cell data. CONCLUSIONS: We propose a G0 arrest transcriptional signature that is linked with therapeutic resistance and can be used to further study and clinically track this state